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Abstract: We study the boundary contribution to the low energy effective action of 
the open string describing the confining flux tube in gauge theories. The form of the 
boundary terms is strongly constrained by the requirement of Lorentz symmetry, which 
is spontaneously broken by the formation of a long confining flux tube in the vacuum. 
\ Writing the boundary action as an expansion in the derivatives of the Nambu-Goldstone 

modes describing the transverse fluctuations of the string, we single out and put in a closed 
form the first few Lorentz invariant boundary terms. We also evaluate the leading deviation 
^ ■ from the Nambu-Goto string produced by the boundary action on the vacuum expectation 

value of the Wilson loop and we test this prediction in the 3d Ising gauge model. Our 
simulation attains a level of precision which is sufficient to test the contribution of this 
term. 
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1. Introduction 

The confinement of quarks is one of the most fundamental features of the strong interactions 
and our understanding of it is far from being complete. However, even if the origin of 
this non-perturbative phenomenon remains unclear, it is still possible to quantitatively 
understand some of its consequences. In particular the vacuum state in presence of a pair 
of quark sources is believed to be dominated, for large quark separations, by the formation 
of a thin flux tube which generates the linear rising of the confining potential. Clearly this 
configuration breaks spontaneously some of the space-time symmetries of the underlying 
gauge theory. The resulting massless Nambu-Goldstone bosons are the natural degrees of 
freedom to be taken into account in the low-energy regime of the theory. They define the 
embedding of the string-like flux tube in the space-time. As we shall see, one may obtain 
important consequences of spontaneous breaking of the space-time symmetries, even if the 
dynamics leading to the formation of the confining flux tube is unknown. 

As first proposed by Liischer, Symanzik and Weisz [1], the effective action of the color 
flux tube connecting a static quark-antiquark pair in a d-dimensional confining Yang-Mills 
theory is a two-dimensional action describing the world-sheet swept out by the flux string 
in its time evolution. The leading term in the infrared limit is a free-field action for the 
d — 2 Nambu-Goldstone modes resulting from the spontaneous breaking of the translation 
invariance in the transverse directions. The ensuing universal string fluctuation effects, 
namely, the leading correction to the linear term in the static potential [2] and the growth 
of the flux tube cross-section with the logarithm of the interquark distance [3] were first 
unambiguously observed many years ago in the Z2 gauge theory in three dimensions [4, 5, 6] 
and more recently in non-abelian Yang-Mills theories [7, 8, 9, 10, 11]. 
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In recent years, the great numerical accuracy of Monte Carlo simulations allowed to 
check and analyze also the subleading corrections of the free-field action. Actually the 
Lagrangian of the low-energy effective theory consists of all terms respecting the internal 
and space-time symmetries of the system. The action of the effective string theory can 
be written as low energy expansion in the number of derivatives of the Nambu-Goldstone 
fields. In particular the first few terms of the effective action S for an open string stretched 
between fixed ends, for instance Polyakov lines, is [12] 

s = s d + ^ J d 2 i [d a x ■ d a x + C2 (d a x ■ d a xf + c 3 (d a x ■ d p xf + ...]+ s b , (1.1) 

where the classical term S c i describes the usual perimeter-area term. The Nambu-Goldstone 
fields Xi(£o, £1) (i = 1, . . . , d — 2) parametrize the displacements orthogonal to the surface 
of minimal area representing the configuration around which we expand and £o>£i are the 
world-sheet coordinates. S b is the boundary contribution characterizing the open string. If 
the boundary is a Polyakov line in the £o direction placed at £i = 0, on which we assume 
Dirichlet boundary conditions -Xj(£o, 0) = 0, the first few terms are 

S b = J d£ [hftX ■ diX + b 2 d 1 d X ■ drfoX + b^X ■ dxXf + ...] . (1.2) 

As first observed in 2004 by Luscher and Weisz [12] the coefficients q and b{ of the above 
expansion should satisfy some consistency constraints; they were obtained by the com- 
parison of the string partition function in different channels ( "open-closed string duality" ) . 
These results were further generalized by Ref. [13]. It was also realized [14, 15, 16] that the 
crucial ingredient of these constraints is the Lorentz symmetry of the underlying Yang-Mills 
theory. Indeed, even if the complete SO(l, d — 1) invariance is broken by the classical con- 
figuration around which we expand, the effective action should still respect this symmetry 
through a non- linear realisation in terms of transverse fields X{ . In this way it was shown 
[15, 16, 17] that the terms with only first derivatives coincide with the Nambu-Goto action 
to all orders in the derivative expansion. The first allowed correction to the Nambu-Goto 
action turns out to be the the six derivative term [13] 

c 4 (d a df)X ■ d a d l3 X^j (d^X ■ d^X) (1.3) 

with arbitrary coefficient C4; however this term is non-trivial only when d > 3. For d = 3 
the first non-trivial deviation of the Nambu-Goto action is an eight-derivative term and 
it has been recently shown [18], using the recursion relations generated by the non-linear 
Lorentz transformations, that it generates a geometric term proportional to the squared 
curvature of the induced metric on the world-sheet. The fact that the first deviations 
from the Nambu-Goto string are of high order, especially in d = 3, explains why in early 
Monte Carlo calculations [19, 20, 21] a good agreement with the Nambu-Goto string was 
observed. The leading deviations of the effective string spectrum were explicitly calculated 
in [22] and they turn out to be consistent with lattice simulations for long closed confining 
strings [23, 24]. 
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In this paper we deal with open effective strings both from the point of view of the- 
oretical investigation and numerical calculations. In [16] it was already observed that in 
the boundary action (1.2) only the 62 term passes the first test of Lorentz invariance (the 
vanishing of b\ was previously proven in [12]). Here we complete the proof by solving 
the recursion relations dictated by the requirement of Lorentz invariance. The resulting 
expansion can be written in a closed form, yielding 



Using the same technique we worked out also the next few boundary invariants. For 
instance, the complete form of that generated by the six derivative term d\d\X ■ d^diX in 



with an arbitrary coefficient b A (this nomenclature will become clear later) . 

The lowest order term in (1.4) represents the first deviation from the Nambu-Goto 
string in d = 3, it is thus important to look for the contribution of this term to physical 
observables in lattice gauge theories. A first test of the effect of this term on the open string 
energy levels has been already performed in [25], in the case of SU(2) gauge theory. In the 
present paper we probe the contribution of this boundary term to the vacuum expectation 
value of Polyakov loop correlators and of Wilson loops, using the three-dimensional Z2 
gauge model. We shall see that this contribution is of the same order of magnitude of 
the subleading effective string corrections to the interquark potential which were recently 
studied both using Polyakov loop correlators [20, 26] and Wilson loops [27] and thus cannot 
be neglected when performing high precision tests of the effective string. We use the data 
obtained from our simulation of Polyakov loop correlators to fit the coefficient 62 of the 
boundary correction. We then use this value to predict the Wilson loop behaviour; we have 
no free parameter left in this comparison, so the remarkable agreement that we find with 
the data represents a strong consistency check of our treatment of the boundary correction 
given in eq. (1.4). 

As it is easy to guess, these boundary corrections are more important in the Wilson 
case than in the Polyakov loop correlators case. In particular, in the Wilson loop case 
they perfectly explain the deviations recently observed in [27] between effective string 
predictions and numerical results. It is possible that they could also explain the subtle, 
higher order, deviations recently observed using Polyakov loop correlators [26], however 
such a test, which would require a much stronger numerical effort is beyond the scope of 
the present paper. 

2. The boundary action 

The confining string action can be thought of as the effective low energy action which is 
obtained integrating over all the massive degrees of freedom of Yang-Mills theory in pres- 




(1.4) 




fiXdfaX (d diX) 2 i (dodiX)Ud 1 xf 
(i + (cU0 2 ) 3 ~ (1 + ^x) 2 ) 4 



(1.5) 



M {l + ^xy? 
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ence of a pair of static sources. As already emphasized in the Introduction, the formation 
of a confining flux tube spontaneously breaks the transverse translational as well as the 
Lorentz (or rotational) invariance of the bulk space-time. This symmetry is not manifest 
in the static gauge, where one uses only the physical degrees of freedom, associated to 
the transverse fluctuations -2Q(£o>£i) (i = 1, ■ ■ ■ d — 2). However the effective string action 
should still respect this symmetry through a non-linear realisation. 



k X i 




Figure 1: Rotation of the Polyakov correlator arount the £o axis . 



To make the discussion concrete and explicit, we specialize now to the case where 
the string action describes the vacuum expectation value of the correlator of a pair of 
Polyakov loops (P(0)*P(R)) parallel to the £o axis and located respectively at £i = 
and £i = R. The string stretched between these two Polyakov loops fluctuates in the 
d — 2 transverse directions, orthogonal to the plane [£o>£i]- Rotate now this system of 
an infinitesimal angle e around the £o axis in the plane as shown in Figure 1. Of 

course this transformation in a Lorentz invariant theory keeps the physics of the system 
invariant; however, it modifies the world-sheet coordinates £o)£i-The linear transformation 
of the transverse fields Xj generated by the rotation should therefore be followed by a 
reparametrization which puts again the system in the static gauge. As a result the string 
action should be invariant under the infinitesimal transformation 

SX j = -eS ij £ 1 -eX i d 1 X j , (i, j = 1, . . . d - 2). (2.1) 

When it is applied to a term with m derivatives and n transverse fields Xj, schematically 
d m X n , it generates terms with the same value of the difference m — n, called "scaling" 
of the given term [18]. When we apply the transformation (2.1) to the boundary term 
dodiX ■ dodiX of (1.2), which has scaling 2, it generates a tower of all possible terms of 
the same scaling value. Thus the Lagrangian density Cb,2 of all boundary terms of scaling 
2 contributing to the boundary action S^p = h% J d£o£b,2 



c b ,2 = Yl [ a kdid x ■ d^xidrX ■ d x x) k + Pkid^oX ■ dxxfidxx ■ d x x) 



fc=0 



(2.2) 



The invariance of £^2 under (2.1) dictates the following recursion relations among the 
and f3k coefficients: 

a k + a k+1 = , (k + l)p k + k/3 k+1 = , a k + f3 k + [3 k+1 = . (2.3) 
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The solution of these recursion relations depends on a single free parameter and yields at 
once Eq.(1.4), which for d = 3 becomes simply 

\2 



d =3 {dpdixy 
(i + {d.xyy 



a, 2 = - ^:r*2 ■ ( 2 - 4 ) 



It is not difficult to enlarge this game to terms of higher scaling. There is no invariant 
term of scaling 3; in fact, it is possible to prove that there are no invariant boundary terms 
of odd scaling. There are two invariant terms of scaling 4. One is the invariant £^4 whose 
form in d = 3 was anticipated in (1.5); we sketch its derivation in appendix B. The other 
is simply (£b,2) 2 , the square of the scaling 2 invariant (2.4). There are five invariants of 
scaling 6, but two of them are simply the product of Lorentz invariants of lower degree. In 
conclusion, the most general boundary action, up to terms of scaling 8, can be written in 
the form 



S h 



3 

b 2 A, 2 + &?> £ M + bf (A )2 ) 2 + £ 6^4,6 + 4 4) (A,2) 3 + 4 5) A,2 A,4 

1=1 



(2.5) 

(i) 

with b k arbitrary coefficients. The general expressions of C\ k for arbitrary space-time 
dimensions d are reported in appendix B. 



3. Effects of the first boundary correction in specific geometries 

In this section we investigate the effect on Polyakov loop correlators and Wilson loop 
expectation values of the first term in the boundary component of the string effective 
action eq. (1.4). We aim for comparison with the lattice simulations discussed in next 
Section, hence we focus on the d = 3 case, where there is a single transverse embedding 
field and the action reduces, as noted in eq. (2.4), to 

^< 2 = 62 / n {d tfJL* = 62 / (Wi*) 2 - 2(d Q d l Xf(d l Xf + ....} (3.1) 

JdT, (1 + (OiX) 2 ) JdT, 

In our simulations we will be able to check the effect of the leading term in the derivative 
expansion of this boundary action, namely 

Sb2 = b ? I (dodiX) 2 (3.2) 

JdT, 

and to fit reliably the value of 62- The first subleading correction in eq. (3.1) appears at the 
same order in the derivative expansion as the leading term arising from the possible scaling- 
weight 4 invariant lagrangian of eq. (1.5); at this order, whose effects in the simulations is 
very tiny, we would therefore have to fit the independent parameter , and this is at the 
moment beyond the reach of our precision. 

Perturbatively, and up to the order discussed above, the correction to the partition 
function for a given world-sheet geometry due to the boundary terms is given by 

^ound^=-^ee(4 1 2 ) ), (3-3) 
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Figure 2: The rectangular Wilson loop geometry that we consider. 



where Zf ree is the partition function obtained using the free action (the first term in eq. 
(1.1)) 

a 



and 



S iree = - JJ(d X) 2 + (d.X) 2 } (3.4) 



b 2 [ lim ****G$& 



(3.5) 



In the last line we used Wick theorem to relate the result to the free propagator (?(£', £) of 
the field X, which depends on the conditions imposed at the boundary <9X. In the following, 
we shall distinguish with a subscript P or W the boundary conditions corresponding to 
Polyakov loop correlators or Wilson loops. 

In [16] the cylinder geometry corresponding to the correlator of two Polyakov loops 
was considered, with the following result: 

<^V - -k^(i^) (3.6) 

where R is the distance between the two Polyakov loops and L the length of the compactified 
time direction (i.e. the inverse temperature). 

Here we consider to be a rectangular Wilson loop of sides R and L, as in figure 2. 
It has four components corresponding to the sides of the rectangle. 

The Green function with Dirichlet boundary conditions on this loop is given by 

r (P t\ 4 A sin sin ggj sin gjfe sin ggl) 

In order to compute eq. (3.5) we focus our attention on one side of the Wilson loop, 
for example the one with coordinates (0, £i), with £i G [0,R\- We have to compute the 
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integral 

" R „ la , a , a Q „ ( fl ^ 4tt 2 ^ m 2 n 2 (cos(^)) 2 



fa rn a 2 ^ 



^— ' / m 2 _j_ n 2 \ 

m ,n=l \ T 2 " ^ ipV 



2tt- ^ nrir 



(3.8) 



i? 2 L 3 



m,n=l ^X 2 " "R 2 " 

With some manipulations we can rewrite the sum in eq. (3.8) as follows: 

v m2 " 2 - 7? 2 y m 2 v 1 1 m2 1 i v^cothT^^ 

m,n=l vl 7 T F' m=l n=l \ \U ^ W j / m=1 

(3.9) 

where we used eq. (A.l) and regularized the divergent sum using eq. (A. 2), taking into 
account that £(0) = —1/2. 

Introducing q = exp(— 2itR/L) we can rewrite the above expression as 

m,n=l IF ^ m=l v ^ 7 I m=l y J 

7T.R 3 Xrn / R' 

(3.10) 

Here we used eq. (A. 2) to regularize the divergent series and eq. (A. 7) to rewrite the 
remaining sum in terms of the Eisenstein series E4. 

Substituting the last expression in eq. (3.8), and using the explicit value £(— 3) = 
1/120, we find thus that the contribution to the boundary term eq. (3.5) from the (0,£i) 
side of the rectangle is given by 



i 



The simmetry of the boundary implies that the contribution from the other side in the 
same direction, the one parametrized by (L,£i), with £1 £ [0,R], must be the same. The 
contribution from the two sides of the rectangle in the orthogonal direction are obtained 
exchanging R and L. Altogether we have thus 

r3 r P / R\ T, f T,\ 1 

(3.12) 



4. Numerical results 

In order to test our prediction for the boundary terms we performed a set of simulations in 
the three-dimensional Z2 gauge model. We used the following strategy. We first evaluated a 
set of Polyakov loop correlators in the low temperature regime. Comparing our results with 
eq. (3.6) we could verify the correctness of the effective string prediction and extract an 
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estimate of the coefficient b 2 . To check the correct scaling behaviour of b 2 we performed this 
analysis for three different values of the bare coupling /3. Then we inserted this estimate for 
62 into eq. (3.12) to predict the boundary correction in the Wilson loop case and compared 
it with the simulation results. In this last comparison there was no free parameter to fit. 
The details of the simulation algorithm that we employed can be found in [28, 27]. As 
a basic update mechanism we used a multispin coded version of the standard metropolis 
algorithm. Let us now describe the various steps of our analysis. 

Polyakov loop correlators 

In order to eliminate the non-universal perimeter and constant terms from the expec- 
tation value of Polyakov loop correlators P(R, L) (where L is the length of the two loops 
and R their distance) we measured the following ratio: 

s P(R+1,L) , . 

Due to the peculiar nature of our algorithm, based on the dual transformation to the 3d 
spin Ising model, this ratio can be evaluated for large values of R and L with very high 
precision. 

The effective string prediction for this observable reads, up to the second loop order, 



ri+ 

'1^2 1! 



R P {R, L) = e~ aL ; ; v 2 *' (1 + F 2 {R + 1, L) + F P {R + 1, L) - F 2 {R, L) - F P (R, L)) , 



(4.2) 

where 77 is Dedekind's function, a is the string tension, F 2 the standard two loop effective 
string contribution to the Polyakov loop correlator (see for instance [20]) 

F 2 (R, L) = ^_(2£ 4 (i —) - El(i—)) (4.3) 

V ' ' 1156i?V V 2R> 2K 2R" V ; 

and Fp = {Sjf^P 1S the leading correction coming from the boundary [16] given in eq. 
(3.6), which we report here for completeness: 

We can use the fact that a is known with very high precision to define a new observable 
which allows to isolate the boundary term: 



v (i L 

R'p(R,L) = e aL m + 2 ' R P (R,L) - F 2 {R+1,L) + F 2 (R,L) (4.5) 

rim) 



From eq. (4.4) we expect for this observable, when L 3> 2R, 

R' p (R,L) = -b 2T ^ (4.6) 

We measured this observable for three values of the bare coupling f3. Table 1 displays the 
values of (3 and the corresponding string tension, as well as the values of L and of the 
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data set 


P 


L 


a 


l/T c 


1 


0.743543 


68 


0.0228068(15) 


5 


2 


0.751805 


100 


0.0105255(11) 


8 


3 


0.754700 


125 


0.0067269(17) 


10 



Table 1: Some informations on the data sample 



R 


data set 1 


R 


data set 2 


R 


data set 3 


11 


0.00611(15) 


16 


0.00515(18) 


20 


0.00446(13) 


12 


0.00431(15) 


17 


0.00395(18) 


21 


0.00355(13) 


13 


0.00290(16) 


18 


0.00334(18) 


22 


0.00304(13) 


14 


0.00182(16) 


19 


0.00207(18) 


23 


0.00226(14) 


15 


0.00159(16) 


20 


0.00132(18) 


24 


0.00190(14) 


16 


0.00092(16) 


21 


0.00119(18) 


25 


0.00162(14) 


17 


0.00088(16) 


22 


0.00133(18) 


26 


0.00127(14) 


18 


0.00062(16) 


23 


0.00083(18) 


27 


0.00134(14) 


19 


0.00020(16) 


24 


0.00071(18) 


28 


0.00104(14) 


20 


0.00016(16) 


25 


0.00088(18) 


29 


0.00074(14) 






26 


0.00067(18) 


30 


0.00060(14) 










31 


0.00078(14) 










32 


0.00062(14) 










33 


0.00053(14) 










34 


0.00033(14) 



Table 2: Numerical results for Rp(R,L)' 



inverse of the critical temperature 1/T C , for the three data sets. We always chose L 3> 1/T C 
to avoid the additional complication of dealing with finite temperature correction terms. 

We report in table 2 the numerical results for the various datasets. We fitted the 
numerical data in table 2 with the function in eq. (4.6), where the only free parameter is 
62 • The results of the fit are given in table 3. For each value of j3 we chose the lowest value 
for R in the fit to be roughly twice the value of the inverse critical temperature. Previous 
analysis on this model [20] showed that this threshold is well within the regime of validity of 
the effective string picture. Looking at table 3 we see that in the fits obtained from all the 
data sets the x 2 values are of order one. This strongly supports the correctness of eq. (3.6), 
obtained in [16]. It is interesting to note that the value of 62 changes dramatically with 
(3. This is due to the fact that 62 has the dimensions of (length) 3 . In fact, a remarkable 
consistency check of our results is that the adimensional combination b2\fo is almost 
constant. This combination can be used to define the continuum limit value of 62 which 
can thus be considered as a new physical scale of the model, on the same ground as the 
string tension a or the glueball mass m g . A slight scaling deviation is still visible between 
the data set 2 and the data set 3 and suggests a conservative estimate of b2\fo ~ 0.032(2) 
for the continuum limit value of hi- It is intersting to notice that our results for 62 in the 
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3d 7L<i gauge model is very similar (but ten times more precise) to the same parameter 
measured in the 3d SU(2) gauge model in [25]. 



data set 


b 2 


62 v 7 ^ 3 


x 2 


1 


7.25(15) 


0.0250(5) 


1.2 


2 


26.8(8) 


0.0289(9) 


1.8 


3 


57.9(12) 


0.0319(7) 


1.3 



Table 3: Some informations on the data sample 



Wilson loops. 

Using the above value of 62 we can now test the leading boundary contribution to the 
Wilson loop reported in eq. (3.12). Since we know a and 62 from independent simulations, 
there are no free parameters left, thus we may perform, so to speak, an "absolute" test of 
our predictions. 

In order to eliminate the non-universal perimeter and costant terms we introduce (as 
in [27]) the following ratio of Wilson loops : 

W(L, R) 



R W {L,R) = 



(4.7) 



W(L + 1,R-1) 

where W(R,L) denotes the expectation value of a Wilson loop of size Rx L (see fig. 2). 
Then in order to eliminate the dominant confining term, using again the fact the a is known 
with very high precision we define 



R W (L, Lu) = R\y(L, uL) — exp{— a(l + L(l — u))} 



(4.8) 



where we introduced the asymmetry parameter u = R/L (with R > L). 

In table 4 we report the numerical results for this observable evaluated for the same 
three values of (3 reported in table 1. 



R data set 1 


R data set 2 


R data set 3 


12 -0.00642(8) 
15 -0.00468(9) 
18 -0.00376(10) 
21 -0.00285(11) 
24 -0.00247(11) 
27 -0.00204(12) 


12 -0.00687(4) 
15 -0.00512(9) 
18 -0.00414(10) 
21 -0.00317(10) 
24 -0.00282(11) 


21 -0.00356(10) 
24 -0.00292(10) 
27 -0.00246(11) 
30 -0.00223(12) 
33 -0.00182(12) 



Table 4: Numerical results for R' W (L, |L) 



In fig.s 3,4 and 5 we plotted the effective string prediction including the one loop and 
two loop corrections from the bulk string action, which we report here for completeness 
(see [27] for a detailed derivation) 



W(L,R) = e -^ L +p(R+ L )+ k 



f](iu) 

~7z 



d-2 
2 



1 + 



C 2 (u) 
aRL 



+ ... 



(4.9) 
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with 



C 2 {u) 



7T \ • 

24 



2(d - 2) u 2 E 4 (iu) 



(d-2)(d-6) 



£ 2 (iu) E 2 (i/n) 



(4.10) 



We also report in the figures the curve obtained adding to the bulk string terms the leading 
correction coming from the boundary, given in eq. (3.12), and the numerical results of our 
simulations. 



1 loop correction 

2 loop correction 
2 loop plus perimeter correction 

numerical data 



22 



24 



26 



Figure 3: R W (L,L±) at /3 = 0.743543 



1 loop correction 

2 loop correction 
2 loop plus perimeter correction 

numerical data 



20 



22 



24 



Figure 4: R' W (L,L^) at (3 = 0.751805 



The agreement between the data and the prediction is impressive and, with the pre- 
cision of our data, the boundary correction is absolutely necessary to correctly describe 
them. Indeed, looking at the similar analysis discussed in [27] (which did not keep into 
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-0.0015 



-0.002 



-0.0025 



-0.003 



-0.0035 



-0.004 



-0.0045 
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1 loop correction 








2 loop correction 








2 loop plus perimeter correction 








numerical data s +■ 



20 22 24 26 28 30 32 34 

Figure 5: R' W (L, §L) at /3 = 0.754700 

account boundary terms) it is clear that the boundary correction discussed in this paper 
exactly fills the gap between data and prediction found, for instance, in fig. 2 of [27]. 



5. Conclusions 



In this paper we used the requirement of Lorentz invariance to constrain the effective 
boundary action of the confining string which describes the behaviour of the Wilson loops 
and Polyakov correlators in any confining gauge theory. We wrote explicitly all the possible 
Lorentz invariants up to terms of scaling 8. We also presented the analytic calculation of the 
contribution of the leading term of the boundary action to the mentioned observables (the 
calculation for the Polyakov correlator was already done in [16]). Although the boundary 
term depends only on (derivatives) of string coordinates evaluated on the boundary, its 
contribution due to quantum fluctuations depends on the whole shape of the observable and 
can be expressed in a closed form in terms of the Eisenstein series E4 . The coefficient 62 in 
front of this boundary term has the dimensions of (length) 3 , thus it defines a new physical 
scale in any confining gauge theory. It is interesting to observe that in the derivative 
expansion of the string action, 62 is the first term where one can find a dependence on the 
gauge group, while all the terms of lower scaling are completely fixed and are the same for 
any 3D gauge theory. We evaluated the coefficient 62 i n the three-dimensional Z2 gauge 
model using the data of Monte Carlo simulations for the Polyakov correlator corresponding 
to three different values of the bare coupling in order to check the correct scaling behaviour 
of 62. We then inserted this estimate to predict the boundary correction in the Wilson loop 
case and compared it with the data of the numerical simulation. Although there are no 
adjustable parameters, the agreement between the prediction and the data is impressive, 
as it is clearly visible in fig.s 3, 4 and 5. 
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A. Useful formulae 

In section 3 we make use of the following series for the hyperbolic cotangent: 

irRL , ,mirR^ L 2 



El ttRL , ,mirR. L z , A . 

— 2 5- = — coth — A -! 
, + it) 2m L 2m 2 

n=l \ L 2 



We also encounter divergent sums which we regularize by means of the zeta function 



oo 



C( a ) = £n-. (A.2) 



n=l 

Dedekind's ^-function is defined as 



V(r)=q^l[(l-q n ) , (A.3) 



n=l 



where (/ = exp(27rir). The Eisenstein series E^fc, with k = 1, 2, . . ., are defined by 



2 



^ 2fe(T) = 1 + C (l-2fc) ^ (J2 ^ l(re)g2re ' (A ' 4) 
where a p (n) is the sum of the p-th powers of the divisors of n: 

a p {n) = Y,mP. (A.5) 

m\n 

In section 3 we encounter the sum 

oo 3 rn oo oo oo 

2 E Si = 2 E E "V" = 2 E ^s(n)5 B , (A.6) 

m=l m=l s=l n=l 

where in the last step we changed summation variables to n = ms and m, so that m is a 
divisor of n, and employed the definition eq. (A.5). From eq. (A.4) we have then 

OO Q m 

ryi ft 

2 E = c(_3) {Ea{t) - l) • (A - 7) 

m=l 

In section 3 the parameter r is in fact iR/L. 
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B. Boundary Lorentz invariants of higher scaling 



In this appendix we list the Lorentz-invariant terms of scaling 4 and 6 which contribute to 
the boundary action of the open string. Since the infinitesimal transformation (2.1) when 
applied to a generic term d m X n keeps invariant the scaling m — n, it suffices to list the 
set of all possible terms of definite scaling with arbitrary coefficients and then solve the 
recursion relations dictated by the request of invariance. For instance in d = 3 there are 
three different sequences of terms of scaling 4, namely 



(B.l) 



£b A [d = 3] = [un(diX) 2n {dldiXf +(3 n (d 1 X) 2n+1 d 2 d 1 X (d d 1 Xy 

n 

+ln (d 1 xf n+2 (dod.xf 

Demanding 5Cb^ = generates the following set of recursion relations 



= (n + 2)a n + (n + l)a n+1 , (B.2) 

= 4a n+ i + (2n + 7)/3„ + (2n + 3)/3„+i , (B.3) 

= Pn+i + (n + 5) 7 „ + (n + 2) 7n+1 . (B.4) 

They can be solved in terms of a single free parameter c = bf , yielding 



a n = (-l)>+l)c,/3 n = -2(-l)"(n+l)(n+2)c, 7 n = ^(-l) n (n+l)(n+2)(n+3)c . (B.5) 



Inserting these solutions in (B.l) gives at once (1.5). It is not difficult to prove that its 
generalization to d space-time dimensions is 



_ d 2 d 1 x • d&hx (ggfrg • d^f + 4 (d^x • d^x) (d^x • &X) 
M ~ i + frx ■ frx (i + ftx • d^f 

(dodiX ■ d!X) 2 [diX ■ d^x + d^x ■ d^x] (d ± x ■ d Q diX) A 



(B.6) 



(1 + diX ■ diX) 6 (1 + diX ■ 5iX) 4 



Using the same procedure it is possible to show that there are in total five Lorentz 
invariants of scaling 6; three of them cannot be written as products of invariants of lower 
scaling. We list them for the sake of completeness: 



14 



c l 6 =7 — o I o ^3 f 3o#l* • ^x \ 9 odiX • a 9ix (i + • ^x) - {d^x ■ d^y 

(1 + d\X ■ OiX) *- L 

- (d^X • diX) (<9 <9iX • d d!X) (dod^ ■ 5iX)(l + diX ■ <9iX) 3 } 

(a^iX • diX) {d Q diX ■ dixf (d^x ■ d G diX) (<9 <9iX ■ d G diX) (d diX ■ diX) 



3 



(1 + diX ■ diXy (1 + diX ■ diX) 

(agftx • fhx) (d d 1 x • d^xf (d 2 d ± x • d^x) {d d 1 x • d^y 

O ■ — q T «J- 



(1 + d x X • diX) 3 (l + diX-3iX) 4 

1 6 {dpdix ■ d x xf 9 (dpdiX • (gggix • gix) (gpdxx • dxX) 2 

(1 + diX ■ diX f ~ ~ ' ' (1 + diX ■ diX) A 

(1 + dxX ■ dxXf + (1 + dxX ■ dxXf 

{d^x ■ d^x) (d d 1 x ■ d!X) 4 
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(i + diX ■ dixy 



(B.7) 



2 _ djd 1 x • gg^x (agAx • ftxy (d^x • ag^x) (giX • d^x) 
bfi l + frx-frx (i + diX-diX) 2 (i + dix-dix) 2 



(<9 3 diX • s diX) (5iX • d 2 d ± x) {dld x x ■ d x x) (dfax ■ ^x) {d^x ■ o^x) 



-6 v u . ; v.. 2 ' + 12 



(1 + diX ■ dxXy (1 + diX ■ dxxy 

| (gfoX • frfrX) (gp^X • ^X) 2 _ (^X • ftX) (do^X • ^X) 5 



(1 + diX ■ diXf (1 + 9iX-9iX) 4 

n {afi&x • d 1 x) 2 (d d 1 x ■ d d 1 x) (d 2 d ± x ■ aft^x) (d d 1 x ■ d^y 

+9 ; ~ z ~t r 9- 



(1 + diX ■ diX) A (1 + diX ■ diX) A 

(dgdiX • a d 1 x) {d 2 d x x ■ d ± x) {d^x • ^x) (d&x • d 1 x) 2 (d d 1 x • a ± xy 



(1 + diX ■ diX) 3 ' (1 + diX • diXf 

(dpdrX • dxX) 2 [(a 2 ^iX • d^X) {dpchX • giX) + (gbgiX • ^X) (gggiX • giX)] 

(i + d x x ■ d x x f 

(9 5iX • diX) 4 [2a 2 5!X • aix + d dix • ao^x] (a aix • diX) 6 



+36^^ ^^-36- 

(1 + diX ■ d x xf (1 + diX ■ 9iX) 6 



(B. 
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6,6 (i + diX ■ dxX) 2 (i + d Y x ■ d!x f 

{dodiX • diX) 4 [dldxX ■ 8 X X + 2d d 1 X • d^X] 

(i + ■ dtX) 5 

+ (1 + • M) 4 

| i {dpd 1 x • c^x) 2 (^giX • d^f _ i (d-fox • d^x) {drfhx ■ d^x) {d^x ■ <hx) 

(1 + d x X ■ diXf (1 + diX • dtX) 3 

(1 + d x X • diX) 4 " (l + d 1 X-d 1 xf 

(B.9) 
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